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ABSTRACT 

The Fermi Gamma-ray Space Telescope recently detected the most energetic gamma-ray burst so far, GRB 
080916C, and reported its detailed temporal properties in an extremely broad spectral range: (i) the time- 
resolved spectra are well described by broken power-law forms over the energy range of 10 keV- 10 GeV, (ii) 
the high-energy emission (at e > 100 MeV) is delayed by w 5 s with respect to the e < 1 MeV emission, and 
(iii) the emission onset times shift towards later times in higher energy bands. We show that this behavior of 
the high-energy emission can be explained by a model in which the prompt emission consists of two compo- 
nents: one is the emission component peaking at e ~ 1 MeV due to the synchrotron-self-Compton radiation of 
electrons accelerated in the internal shock of the jet and the other is the component peaking at e ^ 100 MeV 
due to up-scattering of the photospheric X-ray emission of the expanding cocoon (i.e., the hot bubble produced 
by dissipation of the jet energy inside the progenitor star) off the same electrons in the jet. Based on this model, 
we derive some constraints on the radius of the progenitor star and the total energy and mass of the cocoon of 
this GRB, which may provide information on the structure of the progenitor star and the physical conditions of 
the jet propagating in the star. The up-scattered cocoon emission could be important for other Fermi GRBs as 
well. We discuss some predictions of this model, including a prompt bright optical emission and a soft X-ray 
excess. 

Subject headings: gamma rays: bursts — gamma-rays: theory 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) were only sparsely observed in 
the > 100 MeV r ange, until the Fermi satellite was launched 
on June 1 1 2008 (lAtwood et al.ll2009l) . Now Fermi provides 
extremely broad energy coverage, 8 keV-300 GeV, with high 
sensitivity for GRBs, and is accumulating a wealth of new 
data which open a completely new window on the physics of 
GRBs. The high-energy temporal and spectral data provided 
by Fermi can severely constrain the physical parameters of 
the GRB emission region and the circumburst environment, 
which will lead to a deeper understanding of the central en- 
gine and the GRB progenitors, and will also constrain models 
of high-ene rgy cosmic ray acceleration (for recent reviews, 
iFalcone et a l. 2008; Zhang 2007; Meszaros 2006). 

The GRB detected by Fermi on September 16 2008 
GRB 0809 16C) shows several important new properties 
Abdo et al.ll2009h : 

(i) The time-resolved spectra (with resolution ^ 5 - 50 s) 
are well fitted by an empirical broke n-power-law func- 
tion (the so-called Band function, iBand et aP 119931) 
from 8 keV up to a photon with energy « 13.2 GeV. 

(ii) The e > 100 MeV emission is not detected together 
with the first e < 1 MeV pulse and the onset of the 
e > 100 MeV emission coincides with the rise of the 
second pulse (sa 5 s after the trigger). 

(iii) Most of the emission in the second pulse shifts towards 
later times as higher energies are considered. 
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(iv) The e > 100 MeV emission lasts at least 1400 s, while 
photons with e < 100 MeV are not detected past 200 s. 

(v) The redshift z ~ 4.35 (iGreiner et al.ll2009h and the flu- 
ence « 2.4 x lO"'^ erg cm'^ in the 10 keV - 10 GeV 
range mean that this is the largest reported isotropic 7- 
ray energy release, Zi-^aso — 8.8 x 10^'* erg. 

Some other bursts detected by the LAT detector of Fermi also 
displ ay high-energy la gs, similar to the properties (ii) and/or 
(iii) (lAbdo et al.ll2009h . and then they should be very impor- 
tant to understand the prompt emission mechanism of GRBs. 
We will call the e < 1 MeV emission and the e > 100 MeV 
emission "MeV emission" and "high-energy emission", re- 
spectively. 

A simple physical picture for the property (i) is that the 
prompt emission consists of a single emission component, 
such as synchrotron radiation of electrons accelerated in inter- 
nal shocks of a relativistic jet. In this picture, the peak of the 
MeV pulse could be attributed to the cessation of the emission 
production (i.e., the shock crossing of the shell) and the sub- 
sequent emission could come from the high latitude regions 
of the shell (e.g.. lDermeJl2004l: IZhang et al.ll2006t) . Thus the 
observed high-energy lag for the second pulse (property (iii)) 
requires that the electron energy spectrum should be harder 
systematically in the higher latitude region. This would imply 
that the particle acceleration process should definitely depend 
on the global parameters of the jet, e.g., the angle-dependent 
relative Lorentz factor of the colliding shells, but such a the- 
ory has not been formulated yet. The property (ii) could be 
just due to the fact that the two pulses originate in two internal 
shocks with different physical conditions for which the elec- 
tron energy spectrum of the second intern al shock is harder 
than that of the first one (I Abdo et aO2009l) . 

Another picture is that the prompt emission consists of the 
MeV component and a delayed high-energy component. The 
latter component could be produced by hadronic effects (i.e., 
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photo-pion process and proton synchrotron emission) (e.g., 
lDermeijr2002l : I Asano et alj|2009l) . because the acceleration of 
protons or ions to high energies may be delayed behind the 
electron acceleration. However, it is not clear that hadronic 
effects can reproduce the smooth Band spectrum (see also 
IWan g et al.' 2009). 

In this paper, we discuss a different two-component emis- 
sion picture in which the delayed high-energy component is 
produced by leptonic process (i.e., electron inverse Compton 
scattering) We focus on the effect that the ambient radiation 
up-scattered by the accelerated electrons in the jet can have a 
later peak than that of the synchrotron and synchrotron-self- 
Compton (SSC) emission of the same electrons (correspond- 
ing to the property (iii)) dWang & Meszarosll2006t iFan et alj 
120081: lBeloborodovll2005l) . Provided that the seed photons 
for the Compton scattering come from the region behind the 
electron acceleration region of the jet (see Figure [U, the up- 
scattered high-energy photon field is highly anisotropic in the 
comoving frame of the jet, i.e., the emissivity is much larger 
for the head-on collisions of the electrons and the seed pho- 
tons. As a result, a stronger emission is observed from the 
higher latitude regions, and thus its flux peak lags behind the 
synchrotron and SSC emission. 

Here we propose that the seed photons may be provided 
by the photospheric emission of an expanding cocoon. GRB 
080916C is a long GRB, and it may originate from the col- 
lapse of a massive star. The relativistic jet produced in the 
central region penetrates the star and deposits most of its en- 
ergy output into a thermal bubble, or cocoon, u ntil it breaks 
out of the star (e. g.. Meszaros & Rees 200ll IZhang et"an 
|2003; Miz uta et al.ll2006l) . The cocoon can store an energy 
comparable or larger than the energy of the prompt emis- 
sion of the jet, a nd thus it may make an observable signatur e 
outside the star (Ra mirez-Ruiz et al.ll2002l:lPe'er et al.ll2006 f). 
The cocoon escaping from the star will emit soft X-rays, and 
these can be up-scattered by the accelerated electrons in the 
jet into the high-energy range. The optical thinning of the 
expanding cocoon may be delayed behind the prompt emis- 
sion of the jet, so that the onset of the high-energy emission 
is delayed behind the MeV emission (corresponding to the 
property (ii)). Thus this model has the potential for explain- 
ing the two delay timescales; the delayed onset of the high- 
energy photons (property (ii)) is due to the delayed emission 
of the cocoon, while the high-energy lag within the second 
pulse (property (iii)) is due to the anisotropic inverse Comp- 
ton scattering. We also show that the combination of the 
time-averaged spectra of the SSC and the up-scattered cocoon 
(UC) emission is roughly consistent with the observed smooth 
power-law spectrum (property (i)) (see Figure|2]i. 

As we will explain below, the UC emission is short-lived 
and may not account for the whole high-energy emission, 
which lasts much longer than the MeV emission (property 
(iv)). It is natural that the high-energy emission in later times 
is related to the afterglow, i.e., produced by the external shock 
which propagat es in the circumburst medium. T his possibil- 
ity is studied bv iKumar & Barniol DurarJ (1200 9'). They claim 
that even the onset phase of the high-energy emission is pro- 
duced by the external shock. However, the rise of the flux of 
GRB 0809 16C in the LAT energy range t^) is too steep for 
the external shock to reproduce it. Thus at least the first part 

^ Other ideas for the delayed high-e nergy e mission of GRB 0809 16C 
in the leptonic model are propo sed inlLiI )2009l) . IZou et alj <2009l) . and lFanI 
<2009l) (see also lAoi et al]|2009h . 
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Fig. 1. — Geometry of our model. The jet with opening angle dj ~ 
0.01 and the cocoon with opening angle 9c ~ 0.8 are ejected from the 
collapsing star. The cocoon becomes optically thin at r = Tpj, ~ 2.1 X 

10'* cm £^'52 (Tc/iOy^l^ some time (< 5 s) after the burst trigger (see 
text for detail). The second pulse of the prompt emission is produced by 
the accelerated electrons in the internal shock of the jet at r = r, ~ 2.2 X 
10'* cm r^j (Af,/2 s), and the cocoon X-ray photons are up-scattered by 
the same electrons at r = r, into the GeV range. The jet is assumed to be ob- 
served from the direction of the jet axis, and 9 describes the emitting point of 
the high-latitude emission. The angle x describes the difference of the direc- 
tions of the expansion and the cocoon photon beams, which can be neglected 
for our adopted parameters. 

of the delayed high-energy component of this burst should be 
related to the prompt emission. 

Although some ferm/-LAT GRBs show high-energy lags, 
the observational analysis of the delayed high-energy emis- 
sion is fairly complete only for GRB 0809 16C at the moment, 
and its data is quite extensive. Thus we focus on this burst to 
examine whether the UC emission is viable for its properties 
in the high energy range in this paper However, our modeling 
is generic, and we will show that the UC emission could be 
important for some other Ferm/-LAT GRBs (including short 
GRBs). 

This paper is organized as follows. In § |2l we derive the 
delay time of the onset of the cocoon photospheric emission 
behind the MeV pulse of the jet, and calculate the flux of the 
cocoon emission. In § [3] we constrain the physical condi- 
tions of the emitting region of the jet by assuming that the 
MeV emission is due to the Ist-order SSC radiation, and show 
that the UC flux of the jet and its lag timescale are com- 
patible with the observed GeV emission. In § IH we calcu- 
late the spectrum and lightcurve in a simple numerical man- 
ner and find appropriate values of the physical parameters 
of the cocoon and jet for reproducing the observed spectrum 
at e > 1 MeV and lightcurve of GRB 0809 16C as the com- 
bination of the Ist-order SSC and UC emission. We dis- 
cuss the emission at e < 1 MeV and the thermal emission 
of the jet in § |5] The UC emission for other long and short 
GRBs are discussed in § |6] We summarize our model and 
its implications in § |71 We use keV cm"^ s"' keV~' and 
keV cm"^ s"' for the units of the flux density and flux, respec- 
tively, for comparison with Fermi papers, e.g., Abdo et aP 
(2009). The conversion factors to usual units are useful: 
1 keV cm-2 s"' keV"' = 6.63 x lO'^^ erg cm'^ s"' Hz"' and 
1 keV cm"^ s"' = 1 .60 x 10"^ erg cm"^ s"' . 

2. COCOON DYNAMICS 

The dy namics of the cocoon after e scap ing from the star i s 
studied in lRamirez-Ruiz et alJ (l2002h and lPe'er etaP (Hool). 
Based on these studies, we derive the delay time of the cocoon 
emission behind the prompt emission of the jet and the energy 
flux of the cocoon emission. 

2.1. Time delay of the cocoon emission onset 

The cocoon dynamics is dominated by its radiation energy 
while it is inside the star. We take the total energy and the to- 
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tal mass of the cocoon and the stellar radius as parameters, Ec, 
Mc, and r*, respectively. After the jet breaks out of the star, 
the cocoon expands with the sound speed, Cj = c/\/3, where 
c is the speed of light. The cocoon expands by its own ther- 
mal pressure in the com oving frame as ex pected in th e stan- 
dard fireball model ( Meszaros et al.lll993t iPiran et al.lll993h . 
Its expansion speed in the comoving frame suddenly becomes 
close to c, and then the opening angle of the cocoon mea- 
sured in the central engine frame is given by 9c ~ where 
Fc ex is the Lorentz factor corresponding to Cj. Thus we obtain 
9c ~ y/2j3 ~ 0.8. 'Mors onv etall (120071) have performed de- 
tailed numerical simulations of the cocoon expansion outside 
the star, which show that the opening angle of the expanding 
cocoon is ~ 40°. This is consistent with the above analytical 
estimate. 

The cocoon material accelerates as F cx r and reaches the 
terminal Lorentz factor F^ = Ec/McC^ at 



:r,F, = 5.0x lO'^cmr,.,, ( ^ 



(1) 



Hereafter we adopt the notation Q = 10*2.i cgs units. For 
r > Ks, the cocoon material is dominated by the kinetic energy. 
The radiation stored in the cocoon is released when the opac- 
ity for the electron scattering becomes less than unity. The 
photosphere radius is given by 
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where ct^ is Thomson cross section. 

The cocoon may become optically thin later than the onset 
of the first MeV pulse (i.e., the burst trigger time). Let f = be 
the photon arrival time at the earth emitted at the stellar sur- 
face at the jet breakout. The first MeV pulse produced within 
the jet at r = r,- is observed at 



Af,: = 



2cF- 
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where Tj is the bulk Lorentz factor of the jet and z is the 
source redshift. This timescale is comparable to the angular 
spreading timescale of the pulse, and we can take Af,- ~ 2 s 
for GRB 080916C. The second MeV pulse is observed ~ 5 s 
after the burst trigger, i.e., at f ~ Af, + 5 ~ 7 s. On the other 
hand, the cocoon photospheric emission is observed at 



Afc 



2cF? 
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where we have used z ~ 4.35. The cocoon photospheric emis- 
sion may be observed from t = Ate to f = Ate + At^i, where 
Afrf ^ Ate is the time during which the cocoon will be adi- 
abatically cooled. If internal dissipation occurs in the jet at 
r = ri > Tph making the second MeV pulse within the dura- 
tion of the cocoon emission, the cocoon photons may be up- 
scattered to higher energies by the energetic electrons within 
the dissipation region of the jet, which may be observed along 
with the second MeV pulse. Therefore we require a condition 
Af,- < Afc < Af, -l-5. This condition puts a constraint to the 
physical parameters of the cocoon. 
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We adopt the parameters Ec ~ 10^^ erg and F^ ~ 50 for the 
purposes of calculating the flux of the cocoon emission. 

2.2. Flux of the cocoon emission 

The kinetic energy of the cocoon can be dissipated into ra- 
diation energy, e.g., via internal shocks or magnetic recon- 
nection during the expansion. Since the cocoon material may 
initially have some inhomogeneities caused by the interac- 
tion with the jet and the stellar envelope, on scales ~ ar*, 
where a < 1 is a numerical fac tor, intern al shocks may occur 
at (Rees & Meszaros 1994; Me szaros & Rees 2000) 
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We focus on the case of < < rph, in which the op- 
tical depth for the electron scattering in the dissipation re- 
gion is large and the photons produced by the dissipation un- 
dergo multiple Compton scattering before escape. As a re- 
sult, the emergi ng spectrum is expected to be quasi-thermal 
dPe'er et al. 2006). The condition r, < r j < rph reduces, by 
equations ([T]i, (|2]i, and (|5]) into 



0.01 ( ^ ] <a<OAr:\,. 
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The comoving temperature of the cocoon when its opening 
angle becomes 9c is approximately given by ^ 

1/4 
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where a is the Stefan constant. The comoving temperature 
evolves adiabatically as T' (x r"' for < r < (i.e., in the 
radiation-dominated phase) a nd T' (x r~^/^ at r > r, (i.e., in 
the matter-dominated phase) (iMeszaros et al.l[T993l) . and thus 
the temperature at r = r,; is 



-2/3 



(9) 



If the fraction of the kinetic energy is converted into ther- 
mal energy at r = rj, the corresponding temperature is 



7i' = 



Ec^d/^c 



-,1/4 



27r(l -cosf?c)r^r*Fcfl 
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If the dissipation is not too strong, i.e., e^/ 5, (rd/rsT'^^^, the 
tempe rature at r = is given by T^^ (cf. iRees & MeszarosI 
l2005h . After the dissipation the temperature evolves as T' oc 
r"^/^ since the cocoon shell is still dominated by the kinetic 
energy and the shell spreading effect is negligible up to r ^ 
2r»F^ > Tph. The comoving temperature at the photosphere 
radius is then estimated by 



'ph- 
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The initial temperature of the cocoon may be estimated more ac- 
curately. The temperature of the cocoon ju.st before the jet break-out is 
7j' = [£c/(V,'o)]'/'', where the volume of the cocoon can be approximated 
as y/ ■JT6-^-r^/3 (Bd is the opening angle of the cocoon just before the 
jet break-out.) The cocoon escapes from the star and its comoving volume is 
Vj'^jj Si 27r(l -cosfic)''^ /3 when the opening angle becomes 9^. (The relativis- 
tic contraction effect is negligible.) At that time the temperature is calculated 

as T;; ~ Tl(VrjVl)-^l^ ~ 0.93 T(^^fili\ - To simplify the equations in the 
following, we neglect the weak dependence on 9ci and adopt Tj^,. 
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The photospheric emission of the cocoon is expected to 
have a quasi-thermal spectrum 



(4)' 



(12) 



fore™ <e<e™, 



where SpJJ and are given by 



ell 2.82 ^r;,|^ ^ 1.2 keV r;\f ("^^ ) (13) 



(l+,)3 27r(.p2, 



!31 keV cm-2 s"' keV"^ ^Ist'"*'// > 
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where I'pfJ = e'^'^/h (h is the Planck constant), and we have 
taken the luminosity distance of GRB 0809 16C as t/z. ~ 1 .2 x 
lO^'' cm. These are derived by assuming that the cocoon is 
viewed on-axis, and this assumption is validated in § 13.11 
IPe'er et al.l (2006) carried out detailed numerical calculations 
of the radiative processes in the cocoon and show that the co- 
coon photospheric emission has the above spectral form, with 
f3 -I regardless of the dissipation process if r^/ < rph. The 
cutoff energy e™, is given by ^ 30 x e^^. 

The observation of GRB 0809 16C shows that there is no 
excess from the Band spectrum at the X-ray band, > 10 keV, 
and we obtain a rough upper limit of the cocoon X-ray emis- 
sion £ph^£p° ^ 40 keV cm"^ s"' . This limit leads to another 
constraint on the cocoon parameters. 



r,,„>0.8£2 ( i£ 



(15) 



3. UP-SCATTERED COCOON (UC), SYNCHROTRON, AND SSC 
EMISSION 

We describe here the calculations of the fluxes of the UC, 
synchrotron, and SSC emission of the jet for the second pulse 
of GRB 0809 16C in the internal shock model. The Fermi 
observation allows us to constrain some global parameters of 
the jet (§ 13. It . and we show that the Ist-order SSC radiation of 
electrons accelerated in an internal shock can explain the peak 
energy and the peak flux of the MeV emission (§ 13.21 ). Then 
we show that the same electrons can produce the UC emission 
which lags behind the SSC emission and has the peak energy 
and the peak flux compatible with the observation in the GeV 
range (§[331l. 

3.1. Overall properties of the jet 

We can constrain, from the Fermi observation, the global 
physical parameters of the jet: the bulk Lorentz factor Tj, the 
emission radius of the second pulse r,, and the opening angle 
6j. First of all, from the absence of a 77 absorption cutoff, we 
obtain a lower liniit on the bulk Lorentz factor of the jet, T j > 
870 (lAbdo et al.l2009l) . Since the angular spreading timescale 
of the pulse is Af,- ~ 2 s for the second pulse, similar to the 
first pulse, the emission radius is estimated by 



r; ~ 2cr' 



2 Af,- 
l+z 



2.2 X 10"'cm j 



Af 
2s 



Since GRB 080916C is so bright, it is probable that the jet is 
viewed on-axis, and we adopt this assumption as a simplifica- 
tion (see Figure[T]i. In this case, the cocoon is viewed off-axis, 
since the jet cone will not be filled with the cocoon material. 
The cocoon emission is thus less beamed, but this off-axis 
effect is not significant because the opening angle of the jet 
can be estimated to be small. The isotropic 7-ray energy of 
this burst is 8.8 x lO''"^ erg. To obtain a realistic value of the 
collimation-corrected 7-ray energy, J; 10^' erg, the jet open- 
ing angle is constrained by Oj < 0.015. We adopt dj ~ 0.01, 
and having adopted a nominal value of Tc ~ 50 in accord with 
the observed time delay of the high-energy emission (equation 
|5j, we obtain 



:0.5 < 1. 



(17) 



Thus the off-axis dimming and softening effects are not sig- 
nificant for the cocoon emission. 

3.2. Internal shock model for the MeV emission 

We assume that the jet is dominated by the kinetic en- 
ergy of protons (instead of magnetic energy) and we esti- 
mate the physical parameters of the jet dissipation region 
for the second MeV pulse. In the scenario where the dis- 
sipation is due to internal shocks (Rees & Meszaros 199^ 
(see a lso f anaitescu & Meszaros 2000; Guetta & Granqn 
|200l iGuDta & Zhang. 2007; Kumar & McMahon ,20081 
IBosnjaket al. 2009), the coUisionless shock waves can am- 
plify the magnetic field and accelerate electrons to a power- 
law energy distribution, which then produce synchrotron ra- 
diation and SSC radiation. We find that the Ist-order SSC ra- 
diation can account for the observed MeV emission. At r = r,-, 
the comoving number density of the jet is estimated by 



n = 



LAf,/(l+z)_ Ijl+zf 
Anr^mnC^ 



32nmpC^r^jAtf 



1.7x10^ cm-3L55r-S(^ 



(18) 



where L is the isotropic-equivalent luminosity of the jet, and 
we have used the fact that the dynamical timescale of the emit- 
ting shell measured in the central engine frame is comparable 
to the angular spreading timescale given by Af,/(1 The 
optical depth of the electron scattering is calculated as 



r = ajn'— ~ 2.6 x 10-^ LssFy' 



Af, 
2 s 



(19) 



This indicates that the dissipation region is optically thin, 
which allows the non-thermal synchrotron and SSC emission 
to be observed. 

The internal energy density produced by the internal shock 
is given by u' = n'OpnipC^, where Op is a mean random Lorentz 
factor of the shocked protons and it is of order unity. Assum- 
ing that a fraction eg of the internal energy of the protons is 
carried into the magnetic field, the field strength is estimated 
by 



l/2_3 ( Af/ 



2s 



nll2(J±_\ 



1/2 



B' = (87reBM') ' 2.6 G ^^55 L^; 

(20) 

Assuming that a fraction e^, of the internal energy of the pro- 
tons is given to the electrons, the minimum Lorentz factor of 
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the electrons is given by 

lm = 



p-l irig 



0.4 



(21) 



where p is the index of the electron energy distribution and 
we have assumed p ~ 3.2 (see §|4]for this number). We find 
that the electrons are radiatively cooled significantly within 
the dynamical time for the parameters for fitting the observa- 
tional data (see § |4|i, and the electron energy spectrum aver- 
aged over the dynamical time is expected to be 

N'ij)och~' fo'-7c<7<7., (22) 
\7'^' for7„<7, 

where 7f is the cooling Lorentz factor of the electrons, esti- 
mated by 



7c 



ajB' (l+xi+XiXuc)ri/(cTj) 



(23) 



Here xi and Xuc are the luminosity ratios of the Ist-order SSC 
radiation to the synchrotron radiation and the UC radiation 
to the Ist-order SSC radiation, respectively. The 2nd-order 
SSC radiation luminosity is negligible compared to the Ist- 
order SSC radiation because of the Klein-Nishina suppres- 
sion. Since xi and jCuc are complicated functions of 7^, one 
cannot obtain a simple analytical formula for estimating 7^, 
but we can estimate 7c for the purpose of fitting the data of 
GRB 0809 16C in this model. We find 7,. - 7„, for fitting the 
observational data (see §|4|i, and thus the parameter jci is writ- 
ten as 

4 2 P 
3 p-2 



(24) 



We find that x^c ~ 0.5 in order for the UC emission flux to 
be compatible with the observed flux in the GeV range, and 
thus equation ( l23T l and ( |24] | lead to a formula for estimating 
the cooling Lorentz factor. 
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4x 10-4 
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The synchrotron characteristic energy and the synchrotron 
peak flux (at the synchrotron energy corresponding to 7^) are 
estimated by 



3heB' , 2r 



2 ^-^ J 



4-nmeC '"1+z 



V3e^B'N2Tj(l+z) 



7SC, 



^ 3.4 keV cm- s"' keV l^Tjl j " 9f {^)(i9) 

The synchrotron self-absorption optical depth at £ < £c = 
Smile Hm)^ can be estimated by jMatsumiva & lokal 120031 : 
iToma et al.ll2008b 
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where Ve(x) is the Euler's Gamma function and we have as- 
sumed that 7(. ~ 7,„. Then the ratio of the frequency at which 
= 1 to Ec is given by 

-==^0.20L5, r^,3 0/ \^—^) j . (31) 

The spectrum of the Ist-order SSC emission is described by 

^t^-1/2' 

. (jf) (#) 
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fore < e^*^, 



for e^^ < £ < e^'^ 



for £' 



< £ < £ 



sc 

m 7 



-;V2 



for £„, < £, 



where 



-SC 



(32) 
(33) 



The eF^ peak energy is given by £,^,*^. Equation ( |28] ) and ( |29] l 
imply that the Ist-order SSC radiation can explain the flux of 
the observed MeV emission. 

The 2nd-order SSC emission is significantly suppressed by 
the Klein-Nishina effect. The peak of the 2nd-order SSC 
spectrum is given by 



1 



Ekn = ^jlcmeC~Y^^ — 38 GeV Vj;. 



{—) 
V400/ 



(34) 



3.3. Spectral and temporal behavior of the UC emission 

Here we derive the observed spectrum of the UC emis- 
sion as a function of the polar angle 9 of the emitting re- 
gion on the shell (see Figure [T]). In the comoving frame 
of the jet, the seed cocoon photons are highly anisotropic. 
As a result, the cocoon photons up-scattered by isotropic 
electrons have an anisotropic energy distribution. The spec- 
trum of radiation scattered at an angle 9'^^ relative to the 
direction of the phot on beam in the Thomson s cattering 
regime is given by (Aharonian & Atovan' 1 19811: iBrunettil 
1^01: Wang & Meszaros 2006: Fan et al. 200^ 



.3 X lO^keVcm-2 s"' keV"' L^f F^^^f (l^)(27) 7^4) = ^^rd -cos^) / d^N'ij) dyf,,(l-2y + 2y 



where = [LAf,/(l +z)]/(T jmpC^). The Ist-order SSC char- 
acteristic energy and the SSC peak flux are approximately 



(35) 

where y = e' /[2j^e[{l-cos9[.^)]. This is the scattered radia- 
tion emissivity in the jet comoving frame, A^'(7) is the elec- 
tron energy spectrum, given by equation {22\ . and J'^, is the 
intensity of the seed photons averaged over solid angle, i.e., 
the mean intensity. This equation can be directly derived by 
averaging equation (17) of .Brunetti (.2001.) over the electron 
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directions. When this equation is averaged over cos6>^^, it re- 
duces to the equation for isotropic IC emission (ISari & EsinI 
1200 Ih . The term ^ = 1 -cos^^^, describes the anisotropy of 
the spectrum, and this is due to the fact that the IC scattering 
is strongest for the head-on colHsions between electrons and 
seed photons. This impUes that the UC emission in the ob- 
server frame is stronger from the high-latitude region of the 
shell, so that its flux peak lags behind the onset of the syn- 
chrotron and SSC emission of the same electrons, which have 
isotropic energy distribution in the comoving frame of the jet. 

If the direction of the photon beam is identical to the direc- 
tion of the expansion of the scattering point of the jet, we may 
take 9i,^ = 6' to calculate the observed flux from the scatter- 
ing point with the polar angle 9' with respect to the jet axis 
measured in the jet comoving frame. In this case 



i= 1-COS0^,= 1-COS6'': 



(36) 



In our case, the cocoon is viewed off-axis, and the direction of 
the seed photon beam in the comoving frame is not identical to 
that of the expansion (see Figure [T]!. However, the difference 
is found to be negligible for our adopted parameters. The an- 
gle between the direction of the photon beam and the expan- 
sion direction in the lab frame is x w (fph/rdiOj-O) ~ 10""*, 
since we focus on the region of 6* < Fy' ^ 9j. Our case has 



C = l-cos(6l'±x')' 



(37) 



Since F^x ~ 0.1, and thus the difference between the values 
of ^ for the on-axis and off-axis cases is so small that we can 
use the on-axis equation (|36] |. 

The mean intensity of the cocoon emission measured at r = 
r,- in the central engine frame is calculated as 



1 



4-Trrj l+z 



(38) 



and J'^, = J^JilTj). We perform the integration in equation 
( |35] ) approximately and obtain an analytical form 



where 



;^,(0')~-arae')«X^(e'), 



1 rfi F"" 



(39) 



(40) 



8(e')-- 



for e' < e 



uC 



for e)!^ <£' <e 



uC 



for e"*^ <£'<£, 



fore^^Jj^ < e'. 



uC 



(41) 



Here we define 
and 



(42) 



1+7 1+7 

e^Z' = ic™ ' = ir™ (AT.) 

^ph 2p ph' ^cut 2p cuf 



The eF^ peak energy of the UC emission is given by e„, in 
the case of /? ~ -1 and p ~ 3.2. 

In order to concentrate on the time-averaged spectrum in- 
cluding the high-latitude emission, we calculate the flux of the 
UC emission by neglecting the radial structure of the emitting 
shell for simplicity. In other words, we consider that the emis- 
sion is produced instantaneously in the infinitely thin shell at 
r = r,. This assumption gives us a model lightcurve that re- 
sembles the observed one. Then we obtain a formula for cal- 
culating the flux of the UC emission 



^.(0 = ^ri-;a^(0)^^^|p^, 



where 



e' = (l+z)e 



2F, 



e(t) = Vi 



n \ l+z 



1/2 



(44) 



(45) 



(46) 



and f, is the emission time in the central engine frame (see 
AppendixlAli. 

The peak energy and the peak flux of the UC emission are 
written as functions of the angle parameter q{6) ~T^^f>^- 



pVC _ 2^2 CO €(^) 
^ /m^ph 



160 MeV 



4q 



(1+^)2 

e(0) 



Uoo/ 



IkeV 



(47) 



^580keV cm-^s"' 



ph [1 _|_p25/2]3 

40^2 1 



(i+qr 



-CO pco 
^ph^£ph 



4x 10-4 J V400/ V400/ V30keVcm-2s-i 

where the functions in the brackets [ ] both have values of zero 
at ^ = and q = oo and have peaks of 1 at q = I and ~ 1.4 at 
^ = 2/3, respectively. This means that the UC flux has a peak 
at ^ ~ 1, or ~ Fy', i.e., the peak time of the UC emission 
lags behind that of the SSC emission on the angular spread- 
ing timescale, Af,- ~ 2 s. This is consistent with the observed 
lag of the GeV emission onset behind the MeV emission peak 
of the second pulse of GRB 0809 16C. Here the values of the 
jet parameters r = 4 x lO""* and 7,,, = 7c = 400 are applicable 
for the Ist-order SSC emission of the jet being consistent with 
the observed MeV emission component (see § 13. 2b . This indi- 
cates that the UC emission of the electrons accelerated in the 
internal shock of the jet, emitting the observed MeV emission, 
can naturally explain the observed flux in the GeV range. 

4. CONSTRAINTS ON MODEL PARAMETERS 

In the above two sections we have shown that the up- 
scattered cocoon emission model can explain the observed 
properties (ii) and (iii) of GRB 0809 16C Hsted in §[I] The 
high-energy lag of the second pulse (property (iii)) can be ex- 
plained by the delayed peak time of the UC emission because 
of the anisotropic inverse Compton scattering. If the con- 
straint on the cocoon physical parameters and F^. (equa- 
tion Is]) is satisfied, the cocoon has not released its radiation 
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when the first pulse is produced, and thus the first pulse is 
not accompanied by the UC emission (property (ii)). The ob- 
served high-energy spectrum of the first pulse oc ^-i 63±o.i2 
dAbdo et al. 2009) can be assumed to be produced only by 
the Ist-order SSC emission. This leads to an energy spectral 
index of the accelerated electrons p ^ 3.2 (see equation |32Tl. 
The first and second internal shocks may have similar electron 
acceleration mechanisms, in which case we can take p ^ 3.2 
also for the second MeV pulse. 

In order to confirm whether this model can also explain the 
observed property (i), i.e., the combination of the Ist-order 
SSC and UC emission can reproduce the observed smooth 
power-law spectrum of the second pulse, we perform numer- 
ical calculations of the equation (|44] |. We approximate the 
spectral shape as 



-UC 



1 + 



1 + 



|3+(l/2) 



-3/2 



1 + 



-UC 



(-;j/2)-/i 



(49) 



where we adopt the value s = 2. We can also similarly calcu- 
late the spectrum of the SSC emission (see Appendix|A]i, and 
combine the UC and SSC emission. 

If the flux of the cocoon photospheric X-rays is given, 
i.e., Ec, Tc, r,, and a are given, the fluxes of the UC and 
SSC emission of the jet are determined by the jet parame- 
ters L,r j,Atj,eB, and Cg. Since Af, ~ 2 s is roughly given 
by the observations, and this value is necessary to explain the 
observed high-energy lag timescale, we have four free param- 
eters. On the other hand, we have four characteristic observ- 
ables; the peak fluxes and peak photon energies of the SSC 
component and the UC component. Therefore the jet param- 
eters are expected to be constrained tightly. 

Figure 12] shows the result of the time-averaged spectrum of 
the second pulse for the cocoon parameters 

£,52=1.0, T, = 52, r,,ii=2.5, a = 0.05, (50) 

and /3 = -1.2. These values satisfy the constraints on the co- 
coon parameters, equations (|5]l, (|7]i, and ( fTsl l. The dashed line 
represents the Ist-order SSC component plus the 2nd-order 
SSC component without taking account of Klein-Nishina ef- 
fect, the dot-dashed line represents the cocoon X-ray emis- 
sion, and the dotted line represents the UC component. The 
thick solid line is the combination of all the components tak- 
ing into account the Klein-Nishina effect. The thin solid line 
and the dot-short-dashed lines represent the Band model spec- 
trum fitted to the observed data and the 95% confidence errors, 
respectively (from the Fermi LAT/GBM analysis group, pri- 
vate communication). This figure shows that our model is 
roughly consistent with the observed spectrum at e > 1 MeV. 
Since the observed number flux of the second pulse at ^ 
1 GeV is ~ 1 counts/s, the detection limit for e > 3 GeV is 
roughly given by sF^i^^ ~ lO-' keV cm"^ s"'(e/3 GeV) (The 
effective area of the LAT de tector is appro ximately constant 
for 1 GeV < e < 300 GeV, lAtwo od et a l.1 12009). Thus the 
small bump at ~ 30 GeV cannot be detected. This is con- 
sistent with the non-detection of photons at > 3 GeV. The 
adopted values of the jet parameters are 

L55 = l.l, ry,3 = 0.93, Af, = 2.3s, ei; = 10"^ e<, = 0.4, (51) 



and p = 3.2. The corresponding values of the optical depth 
for electron scattering and the characteristic electron Lorentz 
factors are 



T = 3.5 X 10"^ 7„, = 400, 7r = 390. 



(52) 



Figure |3] shows the results of the multi-band lightcurves for 
the same parameters. Each lightcurve is normalized to a peak 
flux of unity. This clearly displays the lag of the high-energy 
emission onset. We also plot the multi-band lightcurves aver- 
aged in 0.5 s time bins in Figure |4] which corresponds to the 
lightcurv es measured by th e GEM and LAT of Fermi satellite 
(Fig 1 of lAbdo et al.ll2009l) . The observed peak number flux 
of the GeV photons is only ~ 1 counts /bin, so that the rising 
part of the GeV emission at f < 5.5 s could not be detected. 
This is consistent with the onset of the GeV emission being 
delayed behind that of the MeV emission. 

5. DISCUSSION 

5.1. An additional emission component in the sub-MeV 
range ? 

The Ist-order SSC emission and UC emission of the in- 
ternal shock electrons of the jet can explain the main spec- 
tral peak (at e 1 MeV) and the delayed high-energy part 
(at £ > 100 MeV), respectively, of the second pulse of GRB 
080916C (see Figure |2]i. However, we should note a caveat 
for this model. The Ist-order SSC emission has a steep 
low-energy slope, i.e., eF^ oc at e < 1 MeV because of 
strong synchrotron self-absorption effect (see equation [3T]l, 
which is not consistent with the observed low-energy spec- 
ti-um eF, cx 98±o.02 dAbdo et al.ll2009h . On the other hand, 
the 2nd-order SSC emission has the same steep low-energy 
slope, which is consistent with the non-detection of the spec- 
tral bump at e > 1 GeV. If there was no self-aborption effect, 
the 2nd-order SSC emission, even with the Klein-Nish ina ef- 
fect, w ould be so bright as to be detected (see also Wan g et aTl 
120091) . In order to reproduce the observed sub-MeV spectrum, 
we require an additional emission component, with careful 
consideration of the physical conditions of the emitting re- 
gion. An additional sub-MeV emission produced by some 
process near the emission region of the main MeV emission 
of the second pulse would be up-scattered by the electrons 
emitting the MeV emission, like the SSC process without the 
self-absorption, and then the up-scattered flux would easily 
violate the observational upper limits. 

Here we propose, for illustrative purposes, an additional 
sub-MeV emission whose up-scattered emission does not vi- 
olate the observational limits. One can argue that additional 
sub-MeV emission may be produced in another internal shock 
(which we can call a "third internal shock") occurring far in 
front of the second internal shock shell at the same observed 
time. A fraction of the emission from the third internal shock 
propagates backward and is up-scattered by the second shell 
electrons into the high-energy range, but in this case, the elec- 
trons have cooled by the adiabatic expansion and then the up- 
scattered flux is expected to be dimmer than the case men- 
tioned above. 

Figure |5] demonstrates an example of such an additional 
sub-MeV emission. The dashed line is the same as the thick 
solid line in Figure |2] i.e., the time-averaged spectrum of 
the combination of the SSC and UC emission of the second 
internal shock and the cocoon X-ray emission. The dotted 
line represents the synchrotron emission of the third inter- 
nal shock. Here we assume that this third, less energetic 
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internal shock accelerates some fraction / of electron s with 
a power-law energy distr i bution ([B ykov & Meszaros 119961: 
Eichler & WaxmanI 120051: iToma et al.i i2008i: iBosniak et alj 
20091) ' . The adopted parameters are 

L55=0.23, rjj = l-8, Af, = 2.3s, 
eB=1.3xlO-^ e, = 0.35, /= 3x10-1 (53) 

The corresponding values of the scattering optical depth 
for the accelerated electrons and the characteristic electron 
Lorentz factors are 

Tac = 8.2 X lO"'', 7,„ = 1.2 X 10^ 7, = 5.6 x lO"*. (54) 

The SSC emission of the third internal shock electrons are sig- 
nificantly suppressed by the Klein-Nishina effect. The cocoon 
X-rays and the optical synchrotron emission of the second in- 
ternal shock are up-scattered by the third internal shock effi- 
ciently, whose fluxes have been taken into account for estimat- 
ing 7c. The former has the peak flux eF^ ^ 200 keV cm^ s"' at 
e ^ 10 TeV. This may be absorbed by the extragalactic back- 
ground radiation. The latter is shown by the dot-dashed line 
in Figure |5] which may arrive at the earth along with the sec- 
ond pulse in the high-energy range but be so dim as not to be 
detected by the LAT detector. The emission radius of the thkd 
shell is estimated by 

r,„~8.4x lO'^cm. (55) 

Thus the radius where the backward photons from the third 
shell is up-scattered by the second shell is ~ (r, + r,a)/2 ~ 
5.3 x lO'^ cm = ReTir; where /?ex — 2.4. At this radius, the col- 
umn density of the second shell is R^^ times smaller and the 
characteristic electrons Lorentz factors are /?ex times smaller 
because of the adiabatic expansion, so that the optical depth 
for the up-scattering is 7?^^ times smaller than that at r = r,. 
The synchrotron emission from the third shell, up-scattered 
by the expanded second shell, is calculated by using equation 
( [39] l with ^ = 1 - cos(7r - 9'), and the result is shown by the 
thick dashed line in Figure |5] This emission arrives at the 
earth {R^^ - l)Af,- ~ 3.2 s after the onset of the second pulse, 
and may not be detected by the LAT detector. The combined 
spectrum of the second internal shock emission {dashed line), 
the synchrotron emission of the third internal shock (dotted 
line), and the synchrotron emission of the second internal 
shock up-scattered by the third internal shock electrons {dot- 
dashed line) is shown by the thick solid line and it may be 
observed as the second pulse. This is consistent with the ob- 
served spectrum of the second pulse. 

5.2. Internal shock radius vs deceleration radius 

The jet material is finally decelerated by the interaction 
with the circumburst medium, into which the external shock 
propagates. Thus the internal shocks of the jet must oc- 
cur below the deceleration radius of the first internal shock 

' We find it difficult to obtain an additional sub-MeV component by the 
SSC mechanism similar to the main MeV component. A sub-MeV compo- 
nent should be less energetic than the main MeV component, so that it should 
be 7„, < 7c- In this case the flux at ef,^ should make a sub-MeV peak. Then 
e^^ should be about one order of magnitude smaller than and F^^ should be 
comparable with those for the main MeV component, and Ea/sm should be 
smaller than unity. However, such values may not be obtained by equation 
(28), (29) with f|c fSC f^j. < ^^^j gjj ^jjj, _^ _^ 

for any reasonable set of the jet parameter values {L, Fy, A/,-,9p,eB, e^}. 



shell. The first shell is estimated to have the isotropic- 
equivalent kinetic energy E ^ LAf,/(l +z) ~ 4.7 x 10^"* erg 
and the bulk Lorentz factor Tj ^10^ after the production 
of the first MeV pulse. The deceleration radius is then 
rdec ^ [llE/{\6TTT]nmpC^)]^l^ - 10'^ cm n-'/l This would 
be larger than the internal shock radius r,- or ria if n ^ 
1 cm""*. This is in a reasonable range of t ypical values of 
the GRB circumburst m edium density (Pana itescu & Kumad 
I2002t lYost et all 120031) . The deceleration time is esti- 
mated to be fdec ^ ^dec(l +z)/{2cTJ) 9 n"'/^ s. After 
the UC emission ceases (i.e., f > Afc + Af^ ^ 15 s) the 
high-energy emission may be due to the external shock (see 
iKumar & Barniol Duranll2009l) . 

5.3. Photospheric emission of the jet 

We have assumed that the jet is dominate d by kinetic en- 
ergy of baryons instead of magnetic energy. IZhang & Pe'erl 
(2009) have argued that the non-detection of the photo- 
spheric thermal emission of the jet for GRB 080916C may 
exclude a possibility of a baryon-dominated jet and sug- 
gest that a Poynting-flux-dom inated jet is preferred (see also 
iDaigne & MochkovitcE 120021) . They show that if the jet 
is dominated by thermal energy at r = cAf,/(l +z) — 2.8 x 
10^ (Af,/0.5 s) cm, the photospheric thermal emission of the 
jet is so bright that it can be detected, which is inconsistent 
with the observation. However, it is more reasonable that 
the jet has already accelerated up to a Tj w 10^ before reach- 
ing that radius, through the usual adiabatic expansion starting 
from the central region of the progenitor star (r^ ^10^ cm. 
This value corresponds to the gravitational radius of a black 
hole with a mass M = IiMq). Thus the jet is dominated by 
kinetic energy at that radius, and the thermal emission of the 
j et ma y be much weaker than the estimate bv lZhang & Pe'erl 
(l2009l) . 

We assume that the jet is dominated by thermal energy at 
r= re, and it is accelerated by converting the thermal energy 
into kinetic energy. We can derive the flux of the residual 
thermal emission of the jet at the photosphere according to 
the standard fireball model, similar to that of the cocoon (see 
§ lU. Let the initial bulk Lorentz factor of the jet at r = rg be 
> 1 . (It may be possible that a fraction of the black hole 
rotational energy is carried into the kinetic energy of the jet 
by the magnetic field.) Then the bulk Lorentz factor of the 
jet saturates at r^j = r^T j/Tg ~ 10^ cm r^^bT j ^T'"^ . The co- 
moving number density of the jet shell n' = L/{ATir^mpC^T'^) 
and the Thomson optical depth r = arn'r/Tj = 1 define the 
photospheric radius 

^^^-4£^-^-^^^«"'^"^^-^"- ^''^ 

The comoving temperature of the jet shell evolves from Tj ~ 
[L/(47rr^cr^a)]''''* as oc r"' at r^ < r < rgj and as cx r"^/^ at 
r^j < r < Tphj . Then we obtain the peak energy and the peak 
flux of the photospheric emission of the jet shell 

eJ;^ c 43 keV L;l^'^'£r;'/'T%' , (57) 
4f^^460 keV cm-2 s"' L^^Hrf/^T^^l (58) 

(The time-averaged flux is smaller by about a factor of 4.) 
This means that the photospheric thermal emission of the jet 
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is marginally hidden by the prompt MeV emission (see Fig- 
ure |5]l. 

The above analysis does not take account of possible effects 
from the stellar matter The front of the jet is dissipated by the 
interaction with the stellar matter, and the shocked jet matter 
and the shocked stellar matter escape sideways to form the 
cocoon. The shocked jet matter in front of the clean jet at 
the break-out will make an X-ray flash. The brightness of 
the X-ray flash highly depends o n the density profile of the 
star ( Waxman & Meszarosll2003h . We do not discuss this in 
further detail. 

If re-conversion of jet kinetic energy into thermal energy 
due to interaction of the sides of the jet with the high- 
pressure cocoon is strong, the residual thermal emission may 
be brighter than the above estimate. The ratio of thermal en- 
ergy to kinetic energy of the clean, adiabatically expanding jet 
shell at the stellar radius is ~ {r*/rsi)~^^^ ^ 0.05. Thus the re- 
conversion fraction at r < r» is required to be smaller than 5% 
to suppress the thermal flux down to the above estimate. How- 
ever, many numerical simulations with various initial con- 
ditions of the jet and the star show that the re-conversion 
fractions can be much higher than 5% while the jet propa- 
gates in the star (e.g., Zhang et al. 2003; Mizuta et al. 2006t 
iMorsony et al.|2 007: Mizuta & Aloy 2003 and even weU out- 
side the star (iLazzati et al..,200 9). This implies that the pho- 
tospheric thermal emission of the jet for GRB 080916C may 
not be hidden by the SSC emission from the internal shock, 
being inconsistent with the observation. 

It could be possible that the prompt MeV emission of 
GR B 0809 16C is mainly the photospheric ern ission of the 
jet (iMeszaros & Reesll2000HLazzati et al.ll2009h . In this case 
also, the delayed high-energy emission could be due to the 
up-scattering of the cocoon emission off the non-thermal elec- 
trons produced in internal shocks, where the SSC and syn- 
chrotron emission from the internal shocks are dimmer than 
the photospheric emission in the MeV energy range. This sce- 
nario requires different modeling from that in thi s paper . 

Alternatively, as discussed by IZhang & Pe'eii (l2009h . it is 
possible that the jet is initially dominated by magnetic energy 
and the baryons in the jet are accelerated by converting the 
magnetic energy into kinetic energy, so that the jet is domi- 
nated by t he baryon kine tic energy at the internal shock ra- 
dius (e.g., 



iLvubarskv 



Vlahakis & Konigl 20031 iKomissarov e"taLll2009t 
I2009h . The cocoon will be still produced by the 



shocked jet and shocked stellar matter, but a large fraction 
of the cocoon energy will be the magnetic energy. However, 
the magnetic fields could be tangled and dissipated through 
escaping from the jet into the cocoon and interacting with the 
shocked stellar matter and the ambient stellar envelope, so that 
the magnetic energy could become in equipartition to thermal 
(radiation) energy. The dynamics of such a cocoon after es- 
caping from the star is similar to that in the non-magnetized 
jet case. This possibility is speculative, and magnetohydro- 
dynamic simulations of the jet propagation in a star will be 
useful to address it. If this is true, our model based on the 
assumption that the MeV and high-energy emission compo- 
nents are produced by the Ist-order SSC radiation from the 
internal shocks and the up-scattered cocoon radiation, respec- 
tively, will be applicable. It is an open question how bright 
the photospheric thermal emission from the magnetized jet 
can be. 

6. IMPLICATIONS FOR OTHER LONG AND SHORT GRBS 



We have focused on the up-scattered cocoon emission 
for explaining the delayed high-energy emission of GRB 
0809 16C. Since the delayed high-energy emission is observed 
for other bursts, it is useful to discuss whether the up-scattered 
cocoon emission may be important for other bursts. 

Our model involves many parameters; the cocoon parame- 
ters Ec, Fc, r,, and a, and the jet parameters L, 9j, Fy, Af,, 
es, and The cocoon parameters depend significantly on 
the luminosity and the initial opening angl e of the jet and the 
density profile of the progenitor star (e.g.. IMeszaros & ReesI 
2001; Morsony et al. 2007), so that their values for general 
bursts are unclear. Here, for simplicity, we fix the values of the 
cocoon parameters as those adopted above for GRB 0809 16C. 

For further simplicity, we only consider the bursts for which 
the prompt emission in the soft 7-ray band could be produced 
by the Ist-order SSC radiation of the electrons. Indeed, it 
seems difficult to explain the prompt emission of all GRBs in 
the si mple SSC model (Pir an et a l. 2009; Kumar & NarayanI 
120091) . The jet parameters for typical long GRBs are in- 
ferred from various observations; L ~ 10^^ erg s~', Oj :^ 0.1 
(Ghirlandaet al. 2004), and F; ~ 300 (e.g.. lLifliwick & Saril 
l2001tlMolinari e t al. 2007). The other parameters with values 
similar to those adopted for GRB 080916C in this paper, i.e., 
Af,/(1 +z) ~ 0.4 s, €b — IQr^, and — 0.2 may provide 7^. ~ 
7,„ and account for the spectral peak energy e^f ^ 200 keV 
and spectral peak flux ej^pf^ 100 keV cm"^ s"' of typical 
long GRBs with redshifts z ~ 2 in the SSC model (see § I3.2| |. 

In this case we have the value of r^mlc similar to that 
for GRB 080916C, so that the up-scattered cocoon emission 
could be important for those bursts (see equation |48] |. al- 
though we require to study some kinematic effects carefully. 
Now we have 



Fc0,~5> 1, 



(59) 



so that the de-beaming effect of the relativistic motion of the 
cocoon shell for the seed photons of the UC emission is sig- 
nificant. However, this effect may be cancelled out by the 
beaming effect of the relativistic motion of the jet. The inter- 
nal shock radius is r, ~ 2 x 10'^ cm, which is still much larger 
than the photospheric radius of the cocoon rph ~ 2 x 10'"^ cm. 
Then we have the angle between the direction of the photon 
beam and the expansion direction of the scattering point of the 
jet in the lab frame x ~ i^'ph/fdiSj-O) ~ 10"^ (see Figure [U 
and 

r,x~3>l. (60) 

Since Tc9j ^ F^x, the de-beaming and beaming effects may 
compensate each other, so that the seed radiation can have the 
flux and spectrum similar to those for GRB 0809 16C. There- 
fore it is quite possible that the up-scattered cocoon emission 
contributes to typical long GRBs and they account for the de- 
layed high-energy emission. 

The condition T jx > 1 means that x > 6* (since we focus on 
the region of 6* ~ Vf) and S, « 2/(1 +F^6i2). Then these bursts 
will not have the high-energy lags, i.e., the property (iii) (see 
equation |48]|. Since the seed cocoon photons come from the 
front in the jet frame, the UC emission in the observer frame 
is stronger from the line of sight. 

For bursts with smaller Af,/(1 +z) so that r,- < rph, the de- 
lay timescale of the up-scattered cocoon emission behind the 
onset of the internal shock emission should be on the order of 
(1 +z)rphdj/c ^ 10^ s, and thus in this case, the up-scattering 
effect may be important only if the duration of the prompt 
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emission is longer than ~ 10^ s. 

Some short GRB s might originate from the collapses of 
the ma ssive stars (iJaniuk & Progal 120081: IZhang et alj 120091: 
iToma et al. 2005]), and thus they could have the delayed UC 
emission. Even if other short GRBs are produced by the com- 
pact star mergers, it might be p ossible that the jet is a ccompa- 
nied by the delayed disc wind (iMetzger et al.ll2008l) and that 
the emission from the disc wind is up-scattered by the elec- 
trons accelerated in the jet. For either progenitor models, the 
delayed high-energy emission associated with short GRBs, if 
any, would provide an interesting tool to approach their ori- 
gins. 

7. SUMMARY AND FURTHER IMPLICATIONS 

The Fermi satellite provides an unparalleled broad energy 
coverage with good temporal resolution, and it is accumulat- 
ing a steady stream of new data on the high-energy emission 
of GRBs. GRB 080916C, detected by Fermi, showed unex- 
pected properties in the high-energy range, i.e., > 100 MeV, 
which are listed in § [1] These include the new results on 
prompt emission; (i) the time-resolved high-energy spectra 
are described by smooth power-laws up to a photon with en- 
ergy w 13.2 GeV, (ii) the high-energy emission is not de- 
tected with the first MeV pulse (i.e., the first pulse in the 
< 1 MeV range) and begins to be detected along with the 
second MeV pulse, (iii) the second pulse has the later peaks 
in the higher energy bands, and (iv) the high-energy emis- 
sion lasts at least 1400 s, much longer than the duration of 
the MeV emission. We have discussed a model in which 
the prompt emission spectrum consists of an MeV compo- 
nent produced by the SSC emission of electrons accelerated 
in internal shocks in the jet and the high-energy component 
produced by up-scattering of the cocoon X-rays off the same 
electrons, and we have shown that this model can explain the 
above three observed properties (i), (ii), and (iii). The ex- 
panding cocoon may become optically thin some time later 
than the first internal shock of the jet (equation|4|, so that the 
first MeV pulse may not be associated with the up-scattered 
cocoon (UC) emission while the second MeV pulse may be 
associated with it (property (ii), see § |2]l. The UC emission 
has an anisotropic energy distribution in the comoving frame 
of the jet so that the observed UC emission is stronger from 
the higher-latitude region of the shell. This results in the lag 
of the flux peak of the UC emission behind the MeV emission 
onset on the angular spreading timescale (property(iii), see 
§[3j Figure[3] and|4|i. Figure |2] shows that the combination of 
the SSC and UC emission can reproduce the observed high- 
energy spectral data (property (i), see §|4|i. The UC emission 
is short-lived and may not account for the whole high-energy 
emission which lasts longer than the MeV emission (prop- 
erty (iv)). It is natural that the high-energy emission in the 
later times is related to the after glow. This has been shown by 
iKumar & Barniol DuranI (120091) . However, the early portion 
of the high-energy emission should be the UC emission, be- 
cause the external shock cannot reproduce the observed steep 
rise of the flux. 

The observed flux and width of the second pulse of this 
burst indicate that the isotropic-equivalent luminosity of the 
jet L ~ 10^^ erg s"' and the angular spreading timescale of 
the pulse Af, ~ 2 s. The absence of a 77 absorption cutoff 
in the spectra lead s to a jet bulk Lorentz fac tor of Tj ~ 10^ 
dAbdoet all 120091) (see also lAoi et al1l2009l) . (A larger Tj 
would make the internal shock radius larger than the decel- 
eration radius of the shell.) Then in the internal shock model 



with a mean random Lorentz factor of the shocked protons 
9p < 8 (and all of the electrons being accelerated, i.e., / = 1), 
the spectral peak flux at the MeV energy cannot be produced 
by synchrotron emission (see equation |26]l, but can be pro- 
duced by the Ist-order SSC emission of the electrons. The 
observed peak energy and the peak flux of the MeV emis- 
sion constrains the Ist-order SSC emission mechanism requir- 
ing the microphysical parameters as eg — 10"^ <C — 0.4, 
or 7,„ ~ 7c ~ 400. These parameters naturally make the UC 
emission of the same electrons compatible with the observed 
flux at the GeV energy, and Af, ~ 2 s is consistent with the ob- 
served lag of the GeV emission behind the onset of the MeV 
emission. 

In our model, we expect this burst to have had bright syn- 
chrotron emission in the optical band, like the "naked-eye" 
GRB 080319B (Racusin et al. 2008). The Ist-order SSC 
emission peaking at the MeV energy has a steep slope in 
the lower energy range because of a strong synchrotron self 
absorption effect (see Figure |2]l. This effect, along with the 
Klein-Nishina effect, suppresses the 2nd-order SSC emission 
so it is not detected in the £ > 1 GeV range. On the other 
hand, an additional emission component is required to repro- 
duce the observed spectrum in the sub-Me V energy range. We 
propose that another internal shock with the smaller luminos- 
ity L ~ 2 x lO^'* erg s"' and a fraction of electrons accelerated 
/ ~ 3 X lO""' may produce the sub-Me V emission (see § 15.11 
and Figure |5]). This interpretation is somewhat ad hoc but 
plausible, and we argue that it is compatible with the Fermi 
observation. 

Many numerical simulations suggest that the hydrodynam- 
ical jet will store a large fraction of thermal energy because 
of the interaction of the jet with the high-pressure cocoon 
(Morsony et al. 2007; Lazzati et al. 2009). Then the jet will 
emit a bright photospheric thermal radiation, which is not ob- 
served for GRB 080916C (Zhang & Pe'eiil2009.) . It could be 
possible that the prompt MeV emission is mainly the photo- 
spheric emission of the jet (instead of the Ist-order SSC radi- 
ation from the internal shocks), while the high-energy emis- 
sion component is still due to the up-scattered cocoon emis- 
sion from the internal shocks. Alternatively, the jet of GRB 
0809 16C may be initially dominated by magnetic energy. It 
is not clear yet whether a large fraction of thermal energy will 
be stored in the magnetized jet while it propagates in a star. 
The jet could become dominated by baryon kinetic energy at 
large radius, and the cocoon dynamics could be similar to that 
in the non-magnetized jet case, so t hat o ur model in this paper 
would be entirely applicable (see § 15.31 ). 

GRB 0809 16C is the only burst for which the observational 
data is extensive enough and the analysis is fairly complete at 
the moment, so that we have focused on this GRB to discuss 
our model in detail. Since our model involves much informa- 
tion from the observational data and many parameters, it is not 
easy to perform a more general study of the UC emission at 
this stage. However our model is generic, and we have briefly 
shown that it could apply to other Ferm/-LAT GRBs with typ- 
ical parameters L ~ lO^^* erg s~', 9j ~ 0.1, and Tj ~ 300 for 
which the prompt emission in the soft 7-ray range is produced 
by the Ist-order SSC radiation of electrons (see §|6]). We may 
do a comprehensive study of our model when more Fermi- 
LAT GRBs are analyzed. The up-scattering effect could ap- 
ply to short GRBs as well. Investigating the seed photons for 
the delayed high-energy emission of short GRBs would be an 
interesting approach to their origins. 
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A simple prediction of our model is that prompt emission 
spectra of at least some long GRBs would have an excess 
above the Band spectrum around ~ 1 keV due to the cocoon 
photospheric emission, and this excess should have a differ- 
ent temporal behavior from that of the MeV emission. Since 
the Swift XRT detector can observe the GRB prompt emis- 
sion in the 0.3-10 keV range at least 100 s after the triggers, 
the cocoon X-ray emission could in principle be observed if 
its onset times are for, > 100 s. In this case, we could ex- 
amine the up-scattering effect of the cocoon emission if fon 
is smaller than the burst duration T, while the cocoon emis- 
sion may c ontribute to the early X-ray afterglow if fon ^ T 
dPe'er et al.i ,2006). For long GRBs with fon < 100 s, the up- 
scattering effect could be tested by future satellite, e.g., EXIST 
^, which is designed to have a soft X-ray detector and a good 
triggering capability. 

If our model is correct, we can constrain the parameter 
range for which hadronic effects are important on the high- 
energy emission of GRBs, and we can also constrain the mod- 
els of high-energy cosmic ray acceleration. Also, the delay 
time of the onset of the high-energy emission, fon, is directly 
linked to the optical-thinning time of the expanding cocoon, 
which constrains the physical parameters of the progenitor 
star and the cocoon material of GRBs. For GRB 080916C, 
the stellar radius and the total energy Ec and mass Mc of the 
cocoon are constrained to be 0.3 < £'^752(Mc/lO~^M0)^^ < 
0.9 and r,,ii > 0.8 E^. ^^(Mc/IO-^MqT^ (see equations |5] 



and [T5] l. The cocoon energy and the cocoon mass come 
from the jet energy released within the star and the stel- 
lar mass swept by the jet, respectively. These constraints 
therefore provide potential tools for investigating the struc- 
ture of the progenitor star just before the explosion, as well as 
the physical conditions of the jet propagati ng inside the stel- 
lar envelope through either analytical (e .g., iMeszaros & ReesI 
200lt iMatzeri l200l iToma et al 1 120071) or numerical (e.g 



Zhang et al.ll2003t iMizuta et al.ll2006t Morsony et al.ll20(j7h 
approaches. 



We thank K. Murase and the Fermi LAT/GBM group mem- 
bers (especially F. Piron and V. Connaughton, who provided 
the Band model spectrum fitted to the observed data with er- 
rors) for useful discussions. We are grateful to the referee, 
D. Lazzati, for several important comments, and to Y. Z. Fan, 
S. Inoue, K. loka, T. Nakamura, R. Yamazaki, and B. Zhang 
for useful discussions after the submission of this paper We 
acknowledge NASA NNX09AT72G, NASA NNX08AL40G, 
and NSF PHY-0757155 for partial support. XFW was sup- 
ported by the National Natural Science Foundation of China 
(grants 10503012, 10621303, and 10633040), National Basic 
Research Program of China (973 Program 2009CB824800), 
and the Special Foundation for the Authors of National Ex- 
cellent Doctorial Dissertations of P. R. China by Chinese 
Academy of Sicences. 



' |http://exist.gsfc.nasa.gov| 



REFERENCES 



Abdo, A. A., and the Fermi LAT and Fermi GBM Collaborations 2009. 

Science, 323, 1688 
Ahanorian, F. A., & Atoyan, A. M. 1981, Ap&SS, 79, 321 
Aoi, J., Murase, K., Takahashi, K., loka, K., & Nagataki, S. 2009, submitted 

to ApJ(arXiv:0904.4878l 
Asano, K., Inoue, S., & Meszaros, P. 2009, ApJ, 699, 953 
Atwood, W. B., et al. 2009, ApJ, 697, 1071 
Band. D. L., et al. 1993, ApJ, 413, 281 
Beloborodov, A. M. 2005, ApJ, 618, L13 
Bosnjak, Z., Daigne, F, & Dubus, G. 2009, A&A, 498, 677 
Brunetti, G. 2001, Astropart. Phys., 13, 107 
Bykov, A. M., & Meszaros, P 1996, ApJ, 1996, 461, L37 
Daigne, F, & Mochkovitch, R. 2002, MNRAS, 336, 1271 
Dermer. C. D. 2002, ApJ, 2002, 574, 65 
Dermer, C. D. 2004, ApJ, 614, 284 
Eichler, D., & Waxman, E. 2005, ApJ, 627, 861 

Falcone, A. D., et al. 2008, GRB section of the white paper on the status and 

future of ground-based TeV gamma-ray astronomy (arXiv:0810.0520 1 
Fan Y. Z., Piran, T., Narayan, R., & Wei, D. M. 2008, MNRAS, 384, 1483 
Fan, Y. Z. 2009, MNRAS, 397, 1539 

Ghirlanda, G., Ghisellini, G., & Lazzati, D. 2004, ApJ, 616, 331 

Granot, J., Piran, T., & Sari, R. 1999, ApJ, 513, 679 

Greiner, J., et al. 2009, A&A, 498, 89 

Guetta, D., & Granot, J. 2003, ApJ, 585, 885 

Gupta, N., & Zhang, B. 2007, MNRAS, 380, 78 

loka, K., & Nakamura, T. 2001, ApJ, 554, L163 

Janiuk, A., & Proga, D. 2008, ApJ, 675, 519 

Komissarov, S. S., Vlahakis, N., Konigl, A., & Barkov, M. V. 2009, 

MNRAS, 394, 1182 
Kumar, P, & McMahon, E. 2008, ApJ, 384, 33 
Kumar, P, & Narayan, R. 2009, MNRAS, 395, 472 
Kumar, P. & Barniol Duran, R. 2009, arXiv:0905.2417 
Lazzati, D., Morsony, B., & Begelman, M. 2009, ApJ, 700, L47 
Li, Z. 2009, submitted to ApJ I arXiv:0810.2932j 
Lithwick, Y, & Sari, R. 2001, ApJ, 555, 540 
Lyubarsky, Y. 2009, ApJ, 698, 1570 
Matsumiya, M., & loka, K. 2003, ApJ, 595, L25 
Matzner, C. D. 2003, MNRAS, 345, 575 



Meszaros, R, Laguna, R, & Rees, M. J. 1993, ApJ, 415, 181 

Meszaros, R, & Rees, M. J. 2000, ApJ, 530, 292 

Meszaros, R, & Rees, M. J. 2001, ApJ, 556, L37 

Meszaros, R 2006, Rep. Prog. Phys., 69, 2259 

Metzger, B. D., Piro, A. L., & Quataert, E. 2008, MNRAS, 390, 781 

Mizuta, A., Yamazaki, T., Nagataki, S., & Mineshige, S. 2006, ApJ, 651, 960 

Mizuta, A., & Aloy, M. A. 2009, ApJ, 699, 1261 

MoUnari, E., et al. 2007, A&A, 469, L13 

Morsony, B. J., Lazzati, D., & Begelman, M. C. 2007, ApJ, 665, 569 

Panaitescu, A., & Meszaros, P 2000, ApJ, 544, L17 

Panaitescu, A., & Kumar, P 2002, ApJ, 571, 779 

Pe'er, A., Meszaros, P. & Rees, M. J. 2006, ApJ, 652, 482 

Piran, T., Shemi, A., & Narayan, R. 1993, MNRAS, 263, 861 

Piran, T., Sari, R., & Zou, Y. C. 2009, MNRAS, 393, 1 107 

Racusin, J. L., et al. 2008, Nature, 455, 183 

Ramirez-Ruiz, E., Celloti, A.. & Rees, M. J. 2002, MNRAS, 337, 1349 

Rees, M. J., & Meszaros, P 1994, ApJ, 430, L93 

Rees, M. J., & Meszaros, P 2005, ApJ, 628, 847 

Sari, R., & Esin, A. A. 2001, ApJ, 548, 787 

Toma, K., Yamazaki, R., & Nakamura, T. 2005, ApJ, 620, 835 

Toma, K., loka, K., Sakamoto, T., & Nakamura, T. 2007, ApJ, 659, 1420 

Toma, K., loka, K., & Nakamura, T. 2008, ApJ, 673, L123 

Vlahakis, N., & Konigl, A. 2003, ApJ, 596, 1 104 

Wang, X. Y, & Meszaros, P 2006, ApJ, 643, L95 

Wang, X. Y, Li, Z., Dai, Z. G., & Meszaros, P 2009, ApJ, 698, L98 

Waxman, E., & Meszaros, P 2003, ApJ, 584, 390 

Woods, E., & Loeb, A. 1999, ApJ, 523, 187 

Yamazaki, R., loka, K., & Nakamura, T. 2003, ApJ, 593, 941 

Yost, S. A., Harrison, F A., Sari, R., & Frail, D. A. 2003, ApJ, 597, 459 

Zhang, B., et al. 2006, ApJ. 642, 354 

Zhang, B. 2007, Chin. J. Astron. Astrophys., 7, 1 

Zhang, B., & Pe'er, A. 2009, ApJ, 700, L65 

Zhang, B., et al. 2009, ApJ, 703, 1696 

Zhang, W., Woosley, S., & MacFadyen, A. L 2003, ApJ, 586, 356 
Zou, Y. C, Fan, Y. Z., & Piran, T. 2009, MNRAS, 396, 1 163 



12 



Toma et al. 



1 E — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — ' ' ) " 




,8 



8 [keV] 

Fig. 2. — Time-averaged spectrum of the second pulse calculated in the up-scattered cocoon emission model. The Ist-order SSC component plus 2nd-order 
SSC component without including the Klein-Nishina effect (dashed line), the cocoon photospheric emission (dot-dashed line), and the up-scattered cocoon (UC) 
emission (dotted line) are shown. The thick solid line represents the combination of these, taking account of the Klein-Nishina effect, which is roughly consistent, 
at e > 1 MeV, with the Band model spectrum (thin solid line) with 95% confidence errors (dot-short-dashed lines) (from the LAT/GBM group of Fermi). The 
bump at ~ 30 GeV is so dim as not to be detected. The adopted parameters are listed in equation i50\ and j5U . 
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Fig. 3. — Photon number fluxes in several frequency ranges calculated in the up-scattered cocoon emission model. Time resolution is set to be 25 ms. Each 
Ughtcurve is normalized to a peak flux of unity. The peak of the GeV Ughtcurve is delayed behind that of the MeV Ughtcurve. 
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Fig. 4. — Same as Figure|3] but integrated in 0.5 s time bins. The observed number flux of the GeV photons is only ' 
only the delayed peak at5.5s<?<7s could be detected. 
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Fig. 5. — Combined time-averaged emission spectrum of the second and third internal shocks and the cocoon. The dashed line is the same as the thick solid 
line in Figure|2] i.e., the cocoon X-rays + the SSC and up-scattered cocoon (UC) emission of the second internal shock of the jet. A third less-energetic internal 
shock occurring at the same observer time as the second internal shock but far in front of it produces synchrotron emission (dotted line), and this emission is 
up-scattered into the GeV range by the second internal shock shell, arriving at the earth ~ 3.2 s after the onset of the second pulse {thick dashed line). The 
optical synchrotron emission from the second internal shock, up-scattered by the third internal shock electrons, is shown by the dot-dashed line. The combined 
spectrum of the components arriving as the second pulse, i.e., the dashed line + the dotted line + the dot-dashed line, is represented by the thick solid line and it is 
consistent with the Band model spectrum (thin solid line) with 95% confidence errors (dot-short-dashed lines) (from the L AT/GBM group of Fermi). The bump 
at ~ 30 GeV is so dim as not to be detected. The adopted parameters of the third internal shock are listed in equation j53) . 
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APPENDIX 

FLUX OF INSTANTANEOUS EMISSION FROM THIN SHELL 



We can calculate the observed specific flux from the shell expanding toward us relativistically by (iGranot et alJ Il999t 
iWoods & Loebl[T999h 



Fe(f)= ^ / d<j) I d(cose) I drr 



r2(l-/5ycos6l)2' 



(Al) 



where e' = (1 + z)eT f3j cos 9) and t = (1 +z)\t-(r/ c)cos9]. For the instantaneous emission from the infinitely th in shell, the 
comoving emissivity can be approximated by (see also lloka & Nakamurall200j1;lYamazaki et al.|[2003l;lDermerll2004h 



X,(r,0 ^ f^,At'A/S(t-tdS(r-n), 



(A2) 



where At' = rj/cTj and Ar' = ri/Vj are the comoving dynamical timescale and the comoving width of the shell. The delta 
functions represent the approximation of the instantaneous emission at f = in the infinitely thin shell at r = r, . The integration 
can be performed straightforwardly as 

F,(t) = FAtoh ^ 



where 



.3 



Fe(to)= —^»nr-j 



e' = {l+z)e 



l+z. 



IT 



0(t) = V2 



n 



J 
t 

T+~z 



1/2 



For the up-scattered cocoon emission, j'^, is given by equation ( |39] | and this leads to equation (l44l l. 
For the Ist-order SSC emission, the comoving emissivity is given by 

■■T- -fie ), 



47rm(,c2 



where 



/(£')={ 



1/3 

-1/2 

-1/2 



-pIi 



for e' < e^'^' , 



fore^*^' < e' < s^^' , 
fore^*^' < e' < sf^' , 
for ef^' < e', 



and 



This leads to 



sc' _ ljl£ sc sc' . 



1+5 
2r, 



sc SC 



ll£ sc 

2r, " ■ 



Fe(t) = tF,^ 



(A3) 

(A4) 
(A5) 
(A6) 
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